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EXECUTIVE  SUMMARY 


The  Fourier  transform  has  long  been  a  standard  tool  in  signal  analysis.  It  per¬ 
forms  the  fundamental  task  of  analyzing  the  frequencies  which  exist  in  a  given  function 
of  time.  The  Fourier  transform  does  not,  however,  allow  one  to  discern  at  what  time 
particular  frequencies  exist  in  a  signal.  To  overcome  this  predicament,  Gabor  intro¬ 
duced  the  short-time  Fourier  transform  which  permits  a  “time-frequency  analysis”  of 
a  given  signal;  that  is,  it  provides  a  means  of  telling  both  what  frequencies  exist  in  a 
signal  as  well  as  when  they  exist.  This  time-frequency  analysis  led  to  the  development 
of  the  Wigner-Ville  distribution  and  its  connection  with  the  Heisenberg  group. 

More  recently,  the  desire  to  analyze  signals  at  different  levels  of  resolution  has  led 
to  the  notion  of  scale  in  signal  analysis.  While  the  analysis  performed  by  the  Gabor 
transform  is  “fixed  scale,”  scaling  plays  a  fundamental  role  in  the  emerging  field  of 
wavelet  theory.  This  emphasis  on  dilation  naturally  leads  to  the  consideration  of  the 
affine  group  and  the  Mellin  transform. 

In  this  report,  we  begin  with  a  review  of  the  Fourier  transform  and  some  of 
its  characteristics.  We  then  introduce  the  Mellin  transform,  showing  its  relationship 
to  the  Fourier  transform  and  developing  its  analogous  properties.  In  the  following 
section,  the  Paley-Wiener  theorem  is  stated  and  proved,  and  the  Hilbert  transform  is 
introduced.  In  particular,  conditions  which  restrict  the  support  of  a  function’s  Fourier 
transform  are  explored. 

With  this  background,  the  stage  is  set  for  the  introduction  of  the  K-transform 
as  the  composition  of  the  Fourier  and  Mellin  transforms.  We  define  this  transforma¬ 
tion,  derive  its  explicit  expression  as  an  integral  operator,  and  obtain  an  asymptotic 
estimate  for  the  transform  kernel.  Some  properties  of  the  K-transform  are  explored, 
and  the  application  of  this  transform  to  the  work  of  Altes  on  mammalian  hearing  is 
noted. 
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In  the  final  section,  we  develop  a  unified  setting  for  Fourier  and  wavelet  analysis. 
The  connection  of  the  Heisenberg  group  with  the  Gabor  transform  and  Wigner-Ville 
distribution  is  shown,  as  is  the  connection  of  the  affine  group  with  the  wavelet  trans¬ 
form.  We  then  combine  these  notions  in  a  single  setting,  the  affine- Heisenberg  group. 
In  this  setting,  we  close  by  introducing  an  affine  version  of  the  Wigner-Ville  distribu¬ 
tion  in  terms  of  the  K-transform. 
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1.  FOURIER  AND  MELLIN  TRANSFORMS 


For  a  finite  energy  signal,  a  function  /  =  f(x)  in  the  space  L2(Wt)  of  all  square- 
integrable  functions  on  M,  its  Fourier  transform  T(f)  =  /  is  defined  by 

A  f00  •  > 

P7)«)  =  /(f)=/  dx  (U) 

J—oo 

and  its  Fourier  integral  representation  by 

/OO  A 

/«)e2"'f<<{.  (1.2) 

-OO 

Strictly  speaking,  for  T2-functions  these  integrals  must  be  interpreted  as  the  limit 

Ao 

taken  in  norm  (analogously  for  the  inverse  Fourier  transform).  Plancherel’s  theorem 

/oo  A  too 

1/(01  ^  =  /  \f(*)\2dx  (1.3) 

-oo  J  —oo 

for  Fourier  transforms  ensures  that  T  is  an  isometry  from  j L2(R)  onto  L2(®);  on  the 
other  hand,  Parseval’s  theorem 

/OO  too  A  _ 

f(x)g(x)dx  =  /  f(t)g(t)dt  (1-4) 

-oo  J — oo 

for  Fourier  transforms  ensures  that  T  is  a  unitary  mapping  from  L2(lR)  onto  L2(lR). 

The  Fourier  transform  is  a  fundamental  tool  in  signal  analysis.  By  using  a  com¬ 
plex  exponential  as  the  transform  kernel,  this  analysis  allows  the  decomposition  of 
a  signal  into  its  individual  frequency  components  and  establishes  the  relative  inten¬ 
sity  of  each  component. [1]  In  mathematical  discussions  of  the  Fourier  transform,  time 
and  frequency  are  often  distinguished  by  writing  ®  for  the  time  domain  and  ®  for 
the  frequency  domain;  also  we  shall  usually  use  Roman  letters  x,t,  ...  for  the  time 
variable  and  Greek  letters  £,  r/,  ...  for  the  frequency  variable.  Whether  on  the  time  or 
frequency  side,  however,  the  Fourier  transform  exploits  the  additive  property  of  real 
numbers  through  use  of  the  fundamental  property  e2m(x+y)  =  e2vtxe2^ jn  establishing 
the  convolution  property 

?  ’•  f  *9  — *  ? f(x  -  v)g(v)  dv^j  =  (Tf){Tg)  ( 1 .5) 
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as  well  as  the  properties 


^ :/(.  +  «)— *e2™^/  (1.6) 

>(rf)(t  +  TI)  (1.7) 


describing  the  interchange  of  translation  and  modulation.  A  simple  change  of  variables 
yields  the  property 


7 :  /( r.)  — ,  V/M-) 

r  r 

(1.8) 

describing  the  effect  of  dilation.  Two  other  useful  results, 

T  :  f  — *  2wi£J-  f 

(1.9) 

and 

(1.10) 

follow  immediately  from  integration  by  parts.  Lastly, 

r-.r  —  Tf, 

(1.11) 

where  f*(x )  =  /(— x).  In  order  to  see  the  action  of  T  on  some  specific  examples,  let 
us  consider  two  very  important  functions,  the  Gaussian  kernel  and  the  Poisson  kernel. 


Gaussian  kernel 


T :  e 


(1.1 2)(i) 


We  see  this  fact  by  completing  the  square  in  the  exponent  of  the  integrand: 

f°°  e~™2  e~2*ixi  dx  =  r  e-^2+2*^-«2)e-’r<2  dx  =  l°°  e~*{x+ii)2  e~*e  dx 

J  —  oo  J —oo  J —oo 

=  f+°°  e-**e-*e  dy  =  f°°  e~ ^  dy  =  e 

Jii— oo  J— oo 


since 


/°°  _  2 

e  ny  dy  —  1.  The  change  in  the  limits  of  integration  ±  oo  — >  =boo  is 

-oo 

readily  verified  by  applying  Cauchy’s  Theorem  to  the  contour  integral  around  the  box 
in  C  with  vertices  — jR,  R  +  i£,  —R  +  and  then  taking  the  limit  as  R  — >  oo. 


Poisson  kernel 

7  ■  -(— ^)  — >  <r2'l£l  (1.12)(») 

Here  we  observe  that  f(x)  =  1  /tt(1  +  x 2)  is  a  rational  function  whose  denominator  is 
two  degrees  greater  than  the  numerator,  so  we  may  apply  the  Residue  Theorem  for 
rational  functions  (from  complex  analysis)[2]  and  exploit  the  fact  that  /  is  even: 


since  i  is  the  only  pole  of  /  in  C+.  Now  calculate  the  residue 


d  /•  r,  \  {z-iy'W*  e2*Wz  e'2^l 

Res(*i  f(z)e  Kl  )  =  lim -  ■  =  lim . j  -  =  -7—~. 

z-*t  7T(  1  +  Z* )  *-*■*  7r(z  +  l)  Zm 


Therefore, 


/OO 

f(x)e~2*ixi  dx  =  e-2*®. 

-OO 


Let  us  now  introduce  the  Mellin  transform.  Set 

£+  =  max(£,0)  ,  £_=max(-£,0)  (£  €  R)  (1.13) 

and  define  Mellin  transforms 

M+<Kt)  =  r  m(,+u~1/2  ^  M-(j>{t)  =  r  mt-*-1'2  dt  (1.14) 

J— OO  J  —  OO 

Observe  these  are  simply  a  slight  modification  of  the  “classical”  definition  [3] 

f  OO 

m<hz)=  mr'd(,  (us) 

Jo 

where  we  have  set  z  =  2nit  +  1/2  to  make  M+<j)  and  M-<j>  functions  of  t  €  M.  Also 
note  the  obvious  relationship  M+[(t>(— £)]  =  The  corresponding  inverse  Mellin 

transforms  are  given  by 


/OO  - : - - — 

M+m(+lt~1/2  dt  = 

-OO 


/OO  — — 

M.mt- 

-OO 


2tti/— 1/2  j,  _ 


m 

((  >  0), 

0 

(f  <  o); 

m 

(f  <  0), 

0 

((  >  0). 

(1.16) 


(1.17) 


Now  we  point  out  the  fundamental  relationship  between  the  Fourier  and  Mellin 
transforms: 

/OO 

^(e  d( 

-OO 

=  £  me12?"1  j 

(by  mapping  £  — *  log£).  Thus, 

/OO 

*(0#"‘‘~1/a<*e  =  AM(t);  (1.18) 

-oo 

similarly, 

T~X  :  —  f°°  mt^'2  d(  =  (1.19) 

J — OO 


From  this  relationship  we  readily  establish  the  analogue  of  Plancherel’s  theorem  for 
the  Mellin  transforms: 


/oo  roo  rO 

\m?d(=  \m\2d(+  \m\2di 

-oo  JO  J —oo 

/oo  roo 

\eil2<)>(ei)\2  di  +  j  |e€/2^(— e€)| 2 

-OO  J  —  oo 

/oo  roo 

|[e^V(^)]v|2*+  /  |[e<'V(-«{)]'f  * 

-oo  J — oo 

by  Plancherel’s  theorem  for  the  Fourier  transform.  Thus,  from  Eqs.  (1.18)  and  (1.19), 

/oo  roo  roo 

\m\2dt=  \M+m2dt+  /  \M^{t)\2dt.  (1.20) 

-oo  J —oo  J —  oo 

In  the  same  way  we  derive  the  analogue  of  Parseval’s  theorem 

/oo  _  roo  _  roo  _ 

M+mM+m*+  /  M.mM.mdt.  (1.20 

-00  J— 00  J— 00 


In  relating  the  Fourier  transform  to  the  Mellin  transform,  it  is  the  mapping  £  — > 
log£  taking  the  positive  reals  to  M  (and  £  — *  log(— £)  taking  the  negative  reals  to  I1.) 
which  suggests  that  multiplication  on  the  half-line  in  the  Mellin  setting  corresponds 
to  addition  on  the  whole  line  in  the  Fourier  setting,  since  logo;  +  logy  =  log:r?/. 
This  crucial  property  becomes  evident  as  we  work  out  the  following  results  (we  only 
demonstrate  results  for  A4+  here;  those  for  A4_  are  essentially  identical): 


Let 

(0*0X0=  /  0(-)0M  — » 

Jo  U)  LO 

and  observe 

m+  :  0  *  0  — *  m+ f  t°°  0(^-)0M  — )  =  r  r  —  e*H-x/2  # 

\J 0  00  00  )  Jo  Jo  00  LO 

roo  roo  ,  ,  At.) 

=  [  4>(»)f  — 

Jo  Jo  CO 

=  (jf 00  mt2*U-1/2  d()  (jT  Wu)u2irit'1/2  du) . 

Thus 

M+  :  0  *  0  — »  A4+0(f).M+0(f).  (1.22) 


This  establishes  that  the  convolution  property  for  the  Mellin  transform  comes  from 
dilating  one  of  the  convolved  functions  (compare  to  Eq.  (1.5)  for  the  Fourier  transform, 
which  comes  from  translating  one  of  the  convolved  functions). 
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Now  observe  that  for  r  >  0, 

M+:  r-^(r.)  —  /_“  %  =  J"  % 


Consequently, 


/•OO 

, — 2'rrtt  /  ^/^\^27rtt-l/2 

JO 


d{. 


M+  :  r 1/2(/>(r.)  — ►  r-2ir,tA<+<^)  =  e-2^logr)tAf+^(f).  (1.23) 


Also, 

poo 

M+  :  f2”T^(0  =  e2™'(lo^>V(0  — ►  /  ^2ir,^(0e27r,t'1/2  # 

«/o 

poo 

=  /  ^(IK2iri{<+T)-1/2^, 

SO 

Af+:^(0^^+^  +  r).  (1.24) 

Thus,  the  Mellin  transform  sends  dilation  to  modulation,  and  (logarithmic)  modu¬ 
lation  to  translation  (compare  to  Eqs.  (1.6)  and  (1.7),  where  the  Fourier  transform 
sends  translation  to  modulation  and  modulation  to  translation). 

Finally,  for  h  >  0, 

M+  :  ((h~1)l2<j>((h)  — ♦  r  £(*-1>/2^(£fc)£2«*+i/2  a 

Jo  £ 

=  [°°  <l>(u)uj2*itlh<Jl2  ~  =  i  f°°  </)(u>)u)2lTit^h~1^2  dw, 

Jo  hu  h  Jo 

and  so 

*4  :  —  {*>+■!>({)■  (1.25) 

Thus  we  see  that  the  Mellin  transform  acts  upon  exponentiation  as  the  Fourier  trans¬ 
form  acts  upon  dilation  (see  Eq.  (1.8)). 

We  conclude  this  section  with  a  few  more  examples  which  will  be  useful  in  later 
discussion: 

poo  , 

*4 :  rvio  — >  /  r^oe2"'"1^  =  am>(« + a/2«),  (i.26> 

Jo 

poo 

M+  :  e"KI  — ►  /  e-M?**-1'2  d£  =  T(2irit  +  1/2),  (1.27) 

Jo 
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and  from  Bateman  [3], for  Rea  >  0, 

M+  :  e~a(i-K  — ♦  V2a~ST(s)el32/8aD-s[p/V2a\  (s  =  2tt it  +  1/2),  (1.28) 

where  D-s  is  the  parabolic,  cylinder  function 

£>_,(*)  = 

and  'P  is  the  confluent  hypergeometric  series  with  integral  representation 

1  fOO 

*{a,c;x)  =  —  /  e~xtta~x[\  +  t)c~a~l  dt. 

I  [a)  Jo 

We  have  seen  thus  far  a  brief  introduction  to  Fourier  and  Mellin  transforms.  In 
particular,  we  notice  that  T  acts  on  a  function  along  all  of  (—00,00),  while  M+ 
and  A4-  look  only  at  (0,  00)  and  (— 00, 0),  respectively.  In  the  next  section,  we  shall 
characterize  when  exactly  T  takes  a  function  on  (—00,00)  to  (0, 00)  (or  (— 00, 0)),  iti 
order  that  we  might  consider  the  composition  of  these  two  transformations. 
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2.  THE  PALEY-WIENER  THEOREM 


aoo  \!/2 

\f(x  +  iy)\  dx )  <  oo  (i,y£l). 

-oo  / 


Define  the  Hardy  space  H2(C+ )  of  functions  /  as  those  functions  which  are  ana¬ 
lytic  in  the  upper  half-plane  C+  =  {z  €  C  :  Imz  >  0},  and  whose  L 2  norms  along  the 
lines  Imz  =  y  (constant)  are  bounded,  i.e., 

\  1/2 

[/  \f(x  +  iy)\'ds 

y>  0 

Similarly,  define  its  companion  space  H2(C _ )  as  those  functions  which  are  analytic  in 
the  lower  half- plane  C_  =  {z  G  C  :  Imz  <  0},  and 

\  1/2 

vl 

y<0  \J— oo  / 

Functions  /  =  f(x)  which  may  be  analytically  extended  to  C+  or  €_  are  referred  to 
in  signal  analysis  terminology  as  analytic  signals.  The  Paley-Wiener  theorem  charac¬ 
terizes  precisely  when  a  function  belongs  to  //2(C+)  or  H2( C_)  (we  prove  the  case  for 
H2{C+)  here) [4]: 


/  y»oo  \  1/2 

supl  /  |/(x  +  iy)\ dx )  <  oo  (x,yG®). 


(2.1)  Theorem  (Paley-Wiener).  F  G  H2(C+)  if  and  only  if  there  exists  0  €  L2(R)  with 
0  =  0  on  (—oo,0)  such  that 


roo 

F(z)= 

Jo 


(2.1) 


Let  us  precede  the  proof  of  the  theorem  with  the  following  lemma: 


(2.2)  Lemma.  Given  any  zq  in  C+, 


±r 

2  wi  J- 


02tziujx 


oo  (t  0  —  Zq  ) 


du> 


g2TTiZQX 

0 


(X  >  o), 

(x  <  0). 


(2.2) 


Proof,  (a)  Suppose  x  >  0.  Let  F+  =  7+  U  [— i?,  jR],  where  the  curve  7+  is  given  by 
7  +(0)  =  R(el e),  0  <  6  <  7r.  Let  g(u)  =  e2mux.  Then  g  is  analytic  on  €  and 


L[ 

liri  A. 


9(z) 


dz 


.JLf 

2ni  Jo 


2ir i  J- y+  (z  —  Zq) 

1  r  e2iri(Re'elxiRei9 

( Re'0  -  Zq) 


Lf 

lm  A 


JlTtizx 


2i ri  /y+  (z  —  z0 ) 


dz 


de 


< 


1  r * 
2ir  Jo 


e2m  (Rei0)xiRei6 


(Reie  -  Zq) 


d6 
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1  fw  e~21rxRsm&  fl 

| R^e  -  z0\ 
l  r 


<-/ 

2%  Jo 


de  < 


1  r* 
2ir  Jo 


0— 27 rxR  sin# 


R 


=  -! 


—2/KxRsin6 


de 


R/2 

2  W2 


de 


(f  or  R  >  2 1 2T0 1 ) 


=  -/ 


~2/nxR  sin  $ 


de. 


Observe  now  that  sin#  >  20/7T  on  [0,  tt/2],  so 


2  W2 


—2irxR  sin# 


de 


2  W2  o__D/o/i/_\  ...  1  —  e27riJ? 


<l-j 


e-2TrxR(29/-K)  jq  _ 


2ttxR 


as  /l  — *  oo.  By  the  Cauchy  integral  formula,  then 


/*00  g27TtU>!T 


2tt*  J- 


dui 


1  roo  5,(0;) 


=  -/‘ 
27ri  J- 


dio 


(u  —  z0)  27ri  J- 00  (u>  —  z0) 

=  lim  -h (I  9(lLdz+  fR  f^dz)  =  lim  -h  / 

fl-oo  27Ti  VA+  (z  ”  20)  J~K  (2  -  20)  /  H-*°°  Jr 

1  f  9(z) 

2i xi  Jr 


g(*) 

r+  (2  -  z0) 


C?Z 


dz  (yR  >  \z0\) 


/r+  (z  -  z0) 

=  s-(z0)  =  e2*iz°x. 

(b)  Suppose  x  <  0.  Let  T_  =  7_  U  [— R,  R]  (R  >  0),  where  the  curve  7_  is 
given  by  7_(0)  =  Re~ie ,  0  <  0  <  7r.  Again,  let  5(a))  =  e2™"*.  Then  5  is  analytic  on 

C.  Notice  first  that  z0  G  C+  is  always  outside  the  curve  T_  and  g  is  analytic,  so  by 
Cauchy’s  theorem 

9(z) 


L 


T_  (z  -  Zo) 


dz  =  0. 


So, 


1  /  & 

1  *  27r^a: 

1  /  6 

27ri  A-  (2  —  20) 

27T*  /y_  (z  —  z0) 

1  {*  e2vi(Re  ie)x(-iRe~ie)  ^ 


1  r 

2ni  Jo 


{Re~ie  -  z0) 


< 


±f 

2tt  Jo 


1  f*  t2™(Re  ,e)x(-iRe-i<>) 

(. Re- Te  -  z0) 

1  /-7T  p— 2jr|a:|/?sin0  D 

<  —  /  - — — — -de 

2ir  Jo  R/2 

— >  0  as  R 


1  f 7r  g27rxKsin#y^ 


±[ 

2n  Jo 


de<  —  I  T—-e — —7  de 

\Re~%e  -  z0| 


(for  >  2|z0|) 


00 


by  part  (a).  Therefore, 

l  roo  e2iriwx 


/oo  2viwx  1  /  /•  p2lr*z*  /■  e27r«x  \ 

-4 - -du>  =  lim  — — :  (  f  -4 - rdz-  - - r  dz  =  0 

-00  (u>  —  z0)  27t*  v/r_  (z  —  zo)  Ji-  (z  —  zq)  / 
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Hence, 


_1_ 
27 xi 


p2n riujx 

7 - \du  = 

(u  -  Zo) 


^Lnizox 

0 


proving  the  lemma. 


(x  >  0), 

(*  <  o); 


Proof  of  theorem.  Suppose  there  exists  <f>  6  Z/2(®)  with  <f>  =  0  on  (— oo,  0)  such  that 

roc  .  roc 

F(z)  =  I  me2"(‘d(= 

Jo  Jo 

where  z  =  x  -\-  iy  £  C+.  To  see  that  F  is  analytic  on  C+,  we  shall  show  that 

Um 

ZO-+0  Zq 

exists  for  all  z  €  C+ .  So,  choose  any  z  G  C+  and  observe  that 

lim  F{z  +  z«)  -  F(z)  =  Um  J”  d(  -  J”  «>({)e2^ 

zo~+0  Zo  zo~+0  Zo 

We  can  bring  the  limit  inside  the  first  integral  by  the  dominated  convergence  theorem, 
since 


roo  /  r oo  \  1/2  /  i»o o  \  1/2 

I  myi(i-^\d(  <  (/o  woi2#)  (jf 

=  (jf  «OlM)V2(jf°  |e-4'sta(2+»>|^)'/2 

<  oo  for  |Im2o|  <  Imz, 


since  <f>  G  L2(M).  Thus, 

lim  l Xf  +  z°)  ~  F(z) 

2o-*-0  20 


g2iri^o 


*0 


* 


/'OO 

=  /  27T7^(0e2^2df 

Hence,  F  is  analytic  on  C+.  On  the  other  hand,  by  Plancherel’s  theorem, 

/oo  roo  roo 

\F{x  +  iy)\2dx  =  \e~2*wy  <f>(u})\2  dco  <  I  \<f>(u)\ 2  du. 

-oo  JO  Jo 


Thus,  since  <f>  G  L2(m) 

/oo  roo 

|F(:r  +  ?j/)|2 dx  <  /  |<^(w)|2du;  <  oo. 

-oo  «/0 


So,  Fg  #2(C+). 
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Now  suppose  F  £  H2( C+).  Define  Ft,  by 

Fb(z)  =  F(z  +  *6),  b  >  0,  z  (E  C+. 

First  we  wish  to  show  that  Fb  satisfies 

*.,*/- r*V  (2-3) 

2fl7  J-oo  (X  —  2:) 

To  do  this,  we  will  select  a  function  ka  which  approximates  FI,,  invoke  Cauchy’s 
integral  formula  for  ka ,  and  then  show  that  this  carries  over  to  Ft,.  Fix  b  >  0,  and 
define  ka  by 

ka(z)  = -  /  Fb(z  +  x)dx  (2.4) 

a  Jo 

for  a  >  0.  Notice  here  that  since  Fb  is  continuous,  lim  kQ(z )  =  Fb(z).  Let  T+  and  7+ 

ot— ►() 

be  as  in  the  lemma  above.  By  the  Cauchy  integral  formula, 

1  /  \  \  f  ka(u) 

ka{z )  =  — 

2m  Jr j 


I  r+  ( u  —  z) 


du 


for  all  R.  For  any  z  €  C+, 


1 

[  K{u\du 

<Lr 

2ni . 

J' y+  (u  -  z) 

2n  Jo 

ka{Reie)iRei$ 


Consequently, 


<  J_  r  \ka(Reid)\R 
~  2n  Jo  R/2 

1  f  ka  ( u ) 

27T?  />+  (u  —  z) 


du 


(Reie  -  z) 
dO  (f or  R  >  2\z\ 

<  -  r  \UReie)\ d9. 

7 r  Jo 


d6 


Now  observe  that 

p—aRsinO 

\ka(Reie)\  < -  /  \Fb(Reie  +  x)\ dx 

a  Jo 

o-aRsin0  /  ,a  \  1/2  /  ra 


\  1/2 


< 


Jo  dX) 


a  \J 0  /  \Jo 

where  C  =  sup|| Ffe||2  <  00.  Thus 


\  '  r*  p—aRsxnO 

\Fb(Reie  +  x)\2  dxj  <  q1/2  -. 


fc>0 

1  r 


-  r  \ka(Reie)\de  <  -JJ-  ["  e-aRsme  d6  <  -JJJ-  r/2e-aR^d6 
n  Jo  not1/2  Jo  not1/2  Jo 


C 


Ra3/ 2 


(1  —  e  aR ) — >0  as  R — >  00. 
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As  a  result, 


k  (z)  =  —  f  -M.u) 

4  ’  2 iri  7r+  (u  - 


r+  ( u  -  z) 


du 


,im(/  M±du+r  * 

#-+oo\./7,  (u  —  z)  7-fi  (X 


R  kQ(x) 


(X  -  z) 


dx 


-j: 


ka(x) 


dx. 


{x-z) 

To  see  that  this  result  goes  over  to  the  Cauchy  integral  formula  for  Fb,  observe  that 

ka  (  X  ) 


Fb(z)  =  limfcQ(z)  =  lim  f  —  -  dx 
'  a— *0  W  a-0  7-oo  X  -  Z 


-  lim  — - 


[  ka^x)dx  + 

[  ka{x)  dx 

J\x\<n  X  —  Z 

J\x\>n  x  —  z 

(fo7*  all  n) 


=  i  /  mdx+lim.L[  hkidx 

2.TTI  J\x\<n  X  —  Z  a— *0  2tTI  7|a;|>»i  X  —  Z 

since  ka  is  continuous  on  [— n,n].  To  estimate  the  second  term,  note  that, 


f  ka(x) 
J\x\>n  X  — 


dx  <  II ka 


1 


f\x\> 


But  \\Fb\\2  <  oo  for  all  b  >  0.  Thus 


ja;|>n  \x  -  z\2 


\  1/2 

dx  l  < 


const. 


n 


liat 


F,,{z  )  =  lim  (2-.  [  dx  +  lim  2-.  I  bM  ,y 

n->oo  \  2tTI  J\x\<n  X  —  Z  a->0  2l Tl  J\x\>n  X  —  Z 

i  r°°  Jb{x )  j_ 

2iri  j-oo  (i 


{x-z) 

With  this,  the  proof  of  the  theorem  is  nearly  complete.  Observe  that 
v  1  f°°  Fb{u)  f°°  1 

Ft(z)  =  27,LhZT)du  =  L 


-/>»(/: 


027T  iu£ 


2iri{u  —  z) 

by  Lemma  (2.2).  Consequently,  setting  2  =  x  +  ic  €  C+, 


2**^)  (£>  to*-**)* 

UZ  \  /*oo  A 

—  du)dt  =  J  Fb(Oe2^zdt 


0-2*tce2*itx 


f  OO  A 

F6(a;  +  ic)  =  /  Ffc(£)e 
Jo 

Note  now  that  for  all  6,  c  >  0, 

a  Too  foo 

/  n(*)e"2,r4ce"2,r^  dt  =  /  F6(f )c-2"(‘-'cK  dt 

J—oo  J— 00 
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Finally,  set 


/oo  roo  A 

Fb(t  +  ic)t~2^  dt  =  /  Fb+C(t)e~2^  dt  =  F6+c(0- 

-oo  ./ — oo 


*(0 


'e2^Ft(0  (£>0), 

.0  «<0). 

The  previous  result  ensures  that  <f>  is  independent  of  e  >  0.  So,  for  Rez  >  e, 

[°°  <t>tiy**z  dC  =  [°°  F'iQe**'***  d(  =  f°°  Fc(C)e2^{z-lt)  d£ 

Jo  Jo  Jo 

=  FJz  -  it)  =  F(z). 

To  see  that  ^  G  I2(i),  we  note  that 

/oo  roo 

\m\2d(  =  lim  e-**‘W£)|2<J£ 

-oo  o— >-0  Jo 

roo  roo  A 

=  lim/  |e-«V(«|2<if  =  lim/  |fi(£)|2<i£ 

b->0  Jo  b—*Q  Jo 

/oo 

|Ffc(x)|2  dx  (b  y  Plancherel) 

-oo 

<  oo  since  F  G  H2{ C+). 

Hence,  we  have  determined  (j>  as  desired,  and  the  proof  is  complete. 

By  the  same  method  of  proof,  the  following  analogous  result  holds: 


(2.5)  Theorem.  F  G  H2(  C_)  if  and  only  if  there  exists  <j>  G  L2(lR)  with  (f>  —  0  on  (0,  oo) 
such  that 

F(z)  =  f°  m^Zdl  (2-5) 

J  —  oo 


From  Eqs.  (2.1)  and  (2.5)  it  follows  immediately  that 

I2(M)  =  Hl( R)  ©  //!(M),  (2.6) 

where  //±(M)  =  {/  :  /  =  lim  F(x  ±  iy),  3 F  G  H2( C±)},  and  that 

y— +0 

/oo  roo  roo 

|/(x)|!<fa=  /  \h(x)\ 2dx+  \f-(x)\2dx,  (2.7) 

- oo  J —oo  J —oo 

where 

roo  A  rO  A 

f(x)  =  h(x)  +  f4x)=  (2.8) 

JO  J— oo 
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In  the  same  vein,  given  any  function  /  6  L2(m)  we  are  able  to  construct  from  it 
a  function  F  in  H2( C+)  or  H2( C_)  by  means  of  the  Hilbert  transform , 


(«/)(*)  =  PV-  /°°  1^-dy. 

7T  7-00  X  —  W 


-oo  x  —  y 

This  integral  is  a  Cauchy  principal  value,  by  which  we  mean 

pVl  r  fM-dy  =  \im- f  l^Ldy. 

7T  7  — oo  X  J/  e— »0  7r  7|x— j/|>c  it  —  y 

We  proceed  as  follows:  Define  g(  by 


9c(y)  = 


-1 


\y\  > 

elsewhere, 


and  note  that 


r  e-2? riy£ 

lim(7e(£)  =  lim  /  - dy 

'-0  '  e-oA.|>f  V  y 


M>f  y 


=  -2 i  j 

Jo 


.  f°°  sm2irfy 


y 


dy 


__  (  -2i(f ) 


(Z>o) 
2  )  (f  <  0) 


=  -*Vsyn(0  (f  7^  0). 


So, 


(«/)(*)  =  P  V-  ^  dy  =  lim  -  /  iM  dy 

7T  7-00  X  —  y  f->0  7T  7|x-t/|>e  X  —  y 

—  lim  -  f  — — —dy  =  *  gc)(x) 

‘->°TTJ\y\>t  y  e-0  7rV  *  A  ’ 

=  !!5g”(/(0^(0)v(*)  =  |:(/(()(-w)si/»(0)v(^), 


and  thus 


(«/)(*)  =  -<(mvn(())VM- 


We  see  then  by  letting 


/  +  *«/ 


/-*«/ 


that  p,  q  6  L2(ffi)  and 

F+(2)  =  /”  p«)e2"M  =  r 

7  —  oo  7—00  2 

/*oo  A 

=  /  €H2(  €+), 

70 


(2.9) 


(2.10) 


(2.11) 


(2.12) 
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while 


f_(z)  ^  r  m^<a=  r  % 

J —oo  J— oo  Z 

=  f  !((,)<?'*’ df,  e  H2(C-). 

J —oo 

Up  to  this  point,  we  have  concentrated  on  the  L 2  theory  of  the  Fourier  transform, 
utilizing  the  fact  that  T  maps  L2(m)  unitarily  onto  L2(R).  Before  moving  on,  we  make 
some  remarks  about  L 1  theory:  Let  p  and  q  be  conjugate  exponents  ( j;  +  ^  =  1 )  such 
that  1  <  p  <  2.  From  the  Hausdorff-Young  inequality 

ll/ll,  <  ll/ll, 

one  concludes  that 

T  :  LP(R)  — >  Lq{ R). 

For  the  case  p  —  1,  the  stronger  result 

F-.L'iR)— >C0(R) 

holds,  where  C'o(M)  is  the  space  of  continuous  bounded  functions  on  ffi. 

We  now  define  the  space  H1(C+)  in  the  same  fashion  as  we  did  H2(C+):  Let 
Hx( C+)  be  the  space  of  all  functions  F  which  are  analytic  in  the  upper  half-plane  (C+, 
and  whose  L 1  norms  along  the  lines  Imz  =  y  (constant)  are  bounded,  i.e., 

/OO 

\F(x  +  iy)\ dx  <  oo  (x,t/€ffi). 

-OO 

A  function  /  €  /^(M)  is  said  to  be  an  analytic  Ll  -signal  if  /  extends  analytically  to  a 
function  F  E  Hl{ C+).  We  define  the  H1- norm  of  an  analytic  L'-signal  to  be 

/OO 

\F(x  +  iy)\dx, 

-OO 

where  F  is  the  unique  analytic  extension  of  /  to  Hx{ C+).  With  these  definitions  wc 
have  the  following  theorem,  similar  in  content  and  proof  to  the  more  difficult  half  of 
the  Paley-Wiener  theorem: 

(2.13)  Theorem.  If  f  is  an  analytic  L1- signal ,  then 

f  =  0on(— oo,0].  (2.13) 

We  also  state  the  following  results,  due  to  Hardy,  without  proof: 
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(2.M)(z) 


(2.14)  Theorem.  If  f  is  an  analytic  L1- signal ,  then 

jf  1/(01  j<c,II/IIh. 

and 

(jf“l/K)l2y),/2<C’!|l/llH'  (-2.14>(») 

where  C\  and  C2  are  constants  independent  of  f. 


These  results  will  be  necessary  in  our  discussion  of  the  classical  Mellin  transform  and 
convolution  in  the  next  section. 

So,  via  the  Paley- Wiener  theorem  and  the  Hilbert  transform,  we  see  what  is 
necessary  to  restrict  the  support  of  a  function’s  Fourier  transform  to  the  half-line.  The 
following  section  builds  upon  this  notion  as  we  move  on  to  consider  the  composition 
of  the  Fourier  and  Mellin  transforms. 
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3.  THE  K-TRANSFORM 


(3.1)  Definition.  The  K-transform  is  the  composition 

JC  =  (A4+  ©  M-)  o  F  :  L2(R)  — ♦  L2(R)  ©  L2( r)  (3.1) 

o/  the  Fourier  and  Mellin  transforms. 


As  K  is  the  composition  of  unitary  mappings,  it  too  will  be  unitary  and  it  will  take 
L2(m)  onto  L2(R)  ©  L2(m).  There  is  an  explicit  expression  for  1C  as  integral  operator. 
Set 


K(x,  t)  =  (2irix)-^H+1'2'>T(2irit  +1/2).  (3.2) 


(3.3)  Theorem.  On  L2(W)  the  IC-transform  is  given  by 

a  POO 

F+(t)  =  (M+f)(t)  =  /  f(x)K(x,t)  dx,  (3.3)(») 

J  —  OO 

and 

A  roo 

F-(t)  =  =  /  f{x)K(-x,t)dx ;  (3.3)(n) 

J  —  OO 

furthermore,  on  L2(R)  ©  L2(M)  the  inverse  /C-transform  is  given  by 

/OO  _  yOO  _ 

F+(t)K(x,t)  dt  +  /  F-(t)K(—x,t)dt,  (3.3)(m) 

-CO  J  —  OO 

where  in  all  cases  K(x,t )  is  defined  by  (3.2). 


Proof.  In  the  right  half-plane  {z  €  C  :  Rez  >0}  the  T-function  is  defined  by  [5] 

te,soo 

T(z)  =  sz  e~suuz  — ,  (3.4) 

Jo  u 

where  -tt/2  <  Arg(s)  +  6  <  ir/2,  or  A rg(s)  +  6  =  ±ir/2  and  0  <  Rez  <1.  We  will 
consider  only  the  case  where  6  =  0.  Now  by  definition, 


F+(i)  =  (jT 

=  r  /w(  r  ^)  <&. 

J ~oo  \J— oo  £  / 
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The  inner  integral  can  be  evaluated  using  the  definition  given  above  for  the  T-function, 
setting  s  —  2tt ix,  u  =  £,  and  z  =  27rif  +  1/2  (and  S  =  0,  remember).  For  these  choices 


T(2nit  + 


1/2)  =  (2%ix)2*it+1/2  J°°  e~2*ixiZ2*it+l/2  y, 


so 


Thus 


Similarly, 


Thus 


fOO 

K(x,t )  =  (2irix)-^it+1^T(2nit  +  1/2)  =  /  e-2"^w+1/2 

Jo 

/co 

/(a;)(27r*a:)_^27r!t+1/,2^r(27rzt  +  1/2)  da\ 

-OO 

F_(0  -  (M_/)(i)  =  /“  (/"  f(x)e~2™t  dx^™-1'2  d( 

=  /I  /(x)  (  /l  e"W(-02^“1/2  <*r 

=  /“  f{x)  (Jo°°  e2vix^2vit+1/2  y)  dx. 

/OO 

f(x)K(—x,t)dx. 

■OO 


dt 


To  determine  the  inverse  transforms,  observe  that 


Thus 


AT*,  *)  =  [°°  e-2^^2«(+i/2  r  eMxt£-2*it+i/2  J.  =  _*). 

do  Q  Jo  £ 

fOO  A 

AW  =  /  /«)e2"*  it  =  7“'  o  A*;1  o  F+ W 

Jo 

=  J^J00  F+{t)t~2mt-1/2  dtj e2nix^  d£ 

=  J°°  F+(t)(£°  ew^-2wit+1/2  y)  d+ 


Consequently, 


Similarly, 


/OO  1*00  _ 

F+(t)K(—x,—t)dt  =  /  F+(t)K(x,t)dt. 

-OO  J  —  OO 

/-(*)=  /°  f(0^ixid(  =  F-loM-JoF^{x) 

J— OO 

=  /°  (  r  F-(t)C2vit~1/2  dt)e2ntxi  d( 

J — OO  W—  OO  / 
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Consequently, 


=  J°°  eMxi(-0~Mt~1'2  d^j  dt 

=  f  F-{t)[J“  e-2^^-2^+l/2  ^  dt 

/OO  roo  _ 

F-(t)K(x,  - 1 )  dt=  F-(t)K(—x,  t )  dt  . 

-OO  J  —  OO 


This  completes  the  proof. 


Using  the  asymptotic  estimate 

r(27Tif  +  1/2)  ~  v/2^e2^tLo3(2^+l/2)-27r,t-l/2/1  + 


'  '  '  V  12(2ir*<  +  l/2)/ 

for  the  T-function,  we  can  derive  a  corresponding  asymptotic  estimate  for  the  kernel 
K  =  K(x,t): 

K(x,t)  =  (2irix)-^it+1/2)T(2irit  +  1/2) 

\J  2 ffe~ftntt+l/2)l‘Og(2irix)e2ritLog(2nit+l/2)-2irit-l/2  (  j  _j _ \ _ A 


12(2tt*<  +  1/2) 


But  for  z  in  (C, 


and  so 


Log(z)  =  log  \z\  +  i  Arg(z), 


while 


02-FitLog(2irit+\/2)  _  e27rit(log |27r!<+l/2|+iAr5(27rt'«+l/2)) 

=  |2t rit  +  i/2\2vite-2*tAr9(2*it+1/2) 


-(2nit+\/2)Log(2nix)  __  -(27n'i+l/2)  log \2nx\  -(27rii+l/2)t  Arg(2irix) 

C  —  C  c 


l  \  27T«t+l/2  2 


S5n(a?)  e“tV  sgn(x)/4 


Consequently, 


_  /  1  \  27nt+l/2 

ff(aM)  -  ^^V2^wJ  |2ir»i  +  l/2|2ir,<e-,Vsa"^^4 


e7rt(7rs^?i(a;)-2  Arp(27n't+l/2))  /  j  _|_ 


V  12(2™*  +  1/2)  A 

Thus  the  asymptotic  behavior  of  K  depends  critically  on  the  sign  of  both  x  and  t 
because  these  will  determine  the  boundedness  of  the  exponential  term.  First,  suppose 
that  x  >  0.  Then 


Hm  Trt(Trsgn(x)  -  2  Arg(2irit  +  1/2))  =  Hm  7r*(7r  -  2Arg(2'Kit  +  1/2)) 
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7r  —  2  arctan(47rf)  87r/(  1  +  (47rt)2) 

=  lim - - - =  lim - — — - 

*-*  °o  1  /  7T  t  t—>oo  1/7 rt2 


=  lim 


8tt2^2 


<-♦00  (1  4. 167r2^2) 


1/2, 


while 


lim  ■Kt('Ksgn{x)  —  2  Arg(2nit  +  1/2))  =  lim  7r^(7r  —  2  Arg(2nit  +  1/2)) 

t — ►  —  OO  t—*  —  OQ 

=  lim  —  nt(ir  +  2Arg(2wit  +  1/2))  =  —00, 

t — ►OO 

since  Arg(2nit  +  1/2)  >  0  for  all  t  >  0.  Now  suppose  that  x  <  0.  In  the  same  fashion 
as  above,  we  see 

lim  7t£(7t  sgn(x)  —  2Arg(2nit  +  1/2))  =  lim  nt(— n  —  2  Arg(2irit  +  1/2)) 

f — t — ►  OO 


t—+oo 


—  lim  —irtfir  +  2Arg(2nit  +  1/2))  =  —00, 
£—►00 


while 


lim  Trt(7rsgn(x)  —  2  Argi^lnit  +  1/2))  =  lim  nt(— tt  —  2  Arg(2'Kit  +  1/2)) 

=  lim  —nti—n  +  2  Ai'g(2nit  +  1  /2))  =  lim  7ri(7r  —  2  Arg(2nit  +  1  /2))  =  1  /2. 

In  particular,  therefore,  the  term 

^nt(7vsgn(x)— 2  Krg(2'Kit+\/2)) 

in  the  asymptotic  estimate  for  K  is  a  bounded  function  of  x  and  t  and  so  K  itself  is 
uniformly  bounded  in  both  x  and  t. 

As  the  /C-transform  is  less  familiar  than  its  individual  components,  some  examples 
will  perhaps  make  its  fundamental  properties  clearer. 

(3.6)  Examples. 


(i)  The  JC-transform  sends  dilation  to  modulation;  that  is, 

M±  oT  :  a1/2f(ax)  — >  a2*H(M±f)(t)  =  a2*HF±{t). 
This  result  follows  directly  from  Eqs.  (1.8)  and  (1.23),  by  which 

1  i(V 


M+  o  T  :  a1/2f(ax)  — *  M+  ^/(-) 


(3.6)(i) 
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=  a2irit(M+f)(t)  =  aMtF+(t). 


The  case  for  M.~  o  T  is  identical. 


(ii)  For  the  sinc-function, 


M±oT-. 


sin27rx  1 

7 xx  2i xit  +  1/2* 


(3.6)(«) 


This  result  follows  by  direct  calculation.  Letting  \  be  the  characteristic  function  on 
the  interval  [-1,1],  we  obtain 


.  .  _  sin27rx  ,  .  .  ,  , , 

M+of:  —— - >  -M+[x(£)] 


7 rx 


J  ~oo  JO 


^27n‘t+l/2 


2-Kit  +1/2 

Again,  M.-  o  T  follows  identically. 


1 


o  2-Kit  +1/2 


(iii)  For  the  Poisson  kernel, 

M±  o  T  :  —  ("t— r')  — >  T(2wit  +  1/2). 

7T  \1  +  X1  J 

Here  we  reference  Eqs  (1.12)(ii)  and  (1.27),  by  which 

o . T  :  —  X+[e-KI]  =  r(2irt(  +  1/2). 

Once  more,  the  same  follows  for  At-  o  T . 

(iv)  For  the  normalized  Gaussian  with  mean  b  and  variance  a, 


M+  oT: 


ay/2K 

re~2wibt 


(®— 6)2/2o2 


M- 


„  —  27rib£  -i 

(by  Eq.  (1.6)) 


(3.6  ){iii) 


=  M+  — —  JF(e-^/“^)2) \  =  (by  Eq.  (1.8)) 

ay  27T  J 

=  M+[e~  £(***>*«*-*"’«]. 

Now  we  substitute  directly  into  Eq.  (1.28)  with  a  =  |(27ra)2,/?  =  2wib,  and  s  = 
2-Kit  +  1/2  to  obtain 


Note  that  we  may  rewrite  this  result  as 

F+(t)  =  2-M+1/4 )K(-ia,t)<!/(Tcit+l/t,  i;  =£). 

To  determine  A4_  o  .F,  we  utilize  the  identity  M+[<j)(—  £)]  =  so  substituting 

=  —  £  above  yields 

F.(l)  =  =  2-/2(21ra)-r(S)>n|,  i;  =£)  =  F+(i)- 

In  Eq.  (3.6)(i)  we  observe  the  critical  property  that  the  /C-transform  sends  scaling 
to  modulation.  Let’s  investigate  the  action  of  K  upon  translation  (regarding  now  only 
the  case  of  A4+  o  !F): 

M+of:f(v  +  x)^  M+(e2*ixif)  =  [°°  e2^f(C)emt~ 1/2  d(. 

Jo 

Representing  /  by  the  inverse  Mellin  transform  of  F+  (where  F+  ==  ( M+  o  T)f),  this 
expression  becomes 

(M+  o  T)  :  f(v  +  x)  —  jH  ew«  QT  F+(s)r2",s_1/2  ds^3"'*-1/2  df. 

If  Fubini’s  theorem  is  applicable  here,  one  obtains  the  result 

(M+  o  F)  :  /(v  +  x)^f°oo  F+(s )  (^°°  ^ 

Without  further  restrictions  on  F+,  the  inner  integral  is  singular  and  so  only  converges 
for  0  <  Re(2ni(t  —  s))  <  1  (see  Eq.  3.4),  which  is  not  the  case  here  since  s,t  6  ffi. 
But  if,  for  instance,  F+  =  Kf  where  /  is  an  analytic  L1-signal,  then  /( 0)  =  0  and  the 
singularity  of  the  inner  integral  may  be  compensated  for.  If  we  interpret  this  integral 
formally  as  a  gamma  function,  we  then  have  a  kernel  Tx  given  by 

Tx(s)  =  (-2vix)-2™r(2iris) 

such  that 

/OO 

F+(s)Tx(t  —  s)ds  =  TX*  F+.  (3.7) 

-OO 

We  observe  that  the  term  £2’r*(<-s)-1  suggestive  of  the  kernel  of  the  classical  Mellin 
transform 

r  oo 

mm  =  M<K-2vit)  =  /  ^(0ra"t_1  dt  (3.8) 

Jo 
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with  inversion  formula 


(3.9) 


/oo 

- CO 

Notice  that  fXtlfi  exists  if 

Jo  WOIy<°°> 

which  is  precisely  Eq.  (2.14)(i)  when  <j)  =  f  and  /  is  an  //‘-analytic  signal.  This 
classical  Mellin  transform  turns  out  to  be  a  nice  setting  in  which  to  look  at  convolution 
with  respect  to  the  /C-transform.  If  /  and  g  are  //‘-analytic,  then  their  convolution  is 
also  //‘-analytic.  Defining  £  =  DJI  o  we  obtain  the  interesting  result 


consequently, 


&(f  *g)  =m  o  T(f  *g)=  m(fg) 

=  m(J°°  mf(u)t2*iu  du  J°°  mg{v)tMvdv) 
=  m(J°°  f  mf(u)mg(v)^u+v)  dudv') 

=  m(J°°  (y°°  mf{u)mg{v  -  u)  du'jt2™'  dv^j 
=  mom~1  o  (mf*mg)i 

&'■  f  *g  — >  &f  *&g. 


(3.10) 


In  other  words,  .8  sends  convolution  to  convolution.  By  the  same  argument,  one 
deduces  the  similar  result  for  K.: 


W  *  gW  =  (JCf  *  Kg)(t  -  &).  (3.11) 

This  framework  with  convolution  and  the  ^-transform  is  tied  in  directly  with  Altes’ 
pursuit  of  a  signal  representation  which  is  independent  of  time  shift  and  scale  change. [6] 
Let  /  and  g  be  analytic  signals  such  that  f,g  €  L‘( K)  D  L2(R).  Recall  that  f*(x)  = 
f(—x).  By  the  Paley-Wiener  theorem,  if  a  function  h  =  h(x)  is  in  H±( R),  then 
h*  €  //±(K)  since  conjugation  does  not  effect  Th  supported  only  on  (0,  oo)  or  (— oo,  0). 
Observe  that 

*  ■  r —> 

thus 

*:/*—»(*/)*•  (3.12) 
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From  Eq.  (3.10),  then, 


£  :  f  *9 


*/♦(*?)*• 


(3.13) 


Altes  presents  the  — |SUl[2  transform”  (denoted  here  by  0)  to  be  the  magnitude- 
squared  (classical)  Mellin  transform  of  a  signal’s  energy  density  spectrum  [6];  that  is 

This  idea  arises  in  the  context  of  work  to  model  mammalian  hearing.  The  purpose  of 
such  a  transform  is  to  “wipe  out”  time  shift  and  scale  change  in  a  signal,  utilizing  the 
properties  that  \F f\  is  invariant  under  translation  of  /  and  is  invariant  under 
dilation  of  <f>.  We  have  in  fact  captured  this  transform  in  a  more  general  setting  with 
the  ^-transform,  since 

K-77^)  |2  =  |SW  o  T{f  *  g*)\2  =  |*(/  *  g*) |2  =  \af  *  (fig)*]2. 

(Altes’  transform  is  the  special  case  /  =  g.) 

We  note  here  that  Altes,  in  the  same  paper,  presents  one  other  slightly  more 
complicated  transform  involving  composition  of  the  Fourier  and  Mellin  transforms 
(the  “p  —  OT  transform”),  which  preserves  some  phase  information  (unlike  the  |.F|2  — 
|$l|2  transform).  It  may  be  of  interest  to  see  if  this  second  transformation  may  be 
interpreted  in  the  context  of  the  1C  and  A  transforms.  We  leave  this  question  as  an 
open  problem  at  present.  Now  our  interests  turn  to  the  wavelet  transform  and  the 
Wigner-Ville  distribution,  as  we  attempt  in  the  next  section  to  incorporate  these  two 
concepts  with  the  /C-transform. 
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4.  WAVELETS  AND  WIGNER-VILLE 


We  have  so  far  seen  an  introduction  to  the  Fourier  and  Mellin  tranfornis,  some 
basic  properties  of  each,  and  the  presentation  of  the  JC-transform  as  a  composition 
of  these  first  two.  Now  let’s  investigate  how  this  all  might  be  tied  in  with  wavelet 
analysis  and  the  Wigner-Ville  distribution.  First,  recall  that  the  Gabor  transform 
(or  short-time  Fourier  transform)  is  defined  by  the  mapping 

/OO  _ 

/(*)<,(* (4.1) 

-OO 

sending  the  signal  /  into  a  two-dimensional  function  living  in  the  time-frequency 
plane  (t,£).[7]  The  basic  idea  of  the  Gabor  transform  is  that  if  you  want  to  know  the 
particular  time  at  which  particular  frequencies  exist  in  a  signal,  then  look  at  a  small 
piece  of  the  signal  around  the  desired  time  t  and  take  its  Fourier  transform. [8]  Many 
properties  of  the  Fourier  transform  carry  over  to  the  Gabor  transform;  however,  the 
analysis  here  depends  critically  on  the  choice  of  the  window  g  (when  Gabor  introduced 
this  transform,  the  special  case  he  considered  was  a  Gaussian  window). [7]  Given  a 
window  function  g(t),  define  its  bandwidth  A£  by 


llslb 


where  the  square  of  the  denominator  is  the  energy  of  g(t).  Two  sinusoids  will  be 
discriminated  in  frequency  only  if  they  are  more  than  A£  apart. [7]  Thus  A£  is  referred 
to  as  the  frequency  resolution  of  the  short-time  Fourier  transform  analysis  with  window 
g.  Similarly,  define  the  “time  spread”  At  of  g(t)  by 

A,  (r„*2is«i2'*),/2 

IIA 

where  once  again  the  square  of  the  denominator  is  the  energy  of  g(t).  Two  pulses  will 
be  discriminated  in  time  only  if  they  are  more  than  At  apart.  [7]  So,  A*  is  referred  to 
as  the  time  resolution  of  this  same  analysis  with  window  g.  Now,  by  the  Heisenberg 
uncertainty  principle,  the  product  AtA£  of  time  resolution  and  frequency  resolution  is 
bounded  below,  so  one  may  attain  either  arbitrarily  good  time  or  frequency  resolution, 
but  not  both.  Once  the  selection  of  g  has  been  made,  the  resolution  over  all  time  and 
frequency  is  fixed.  Therefore,  an  alternative  to  the  Gabor  transform  is  necessary  if  we 
desire  varying  levels  of  resolution  in  the  analysis  of  a  given  signal.  Recently,  wavelet 
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analysis  has  been  presented  as  one  such  alternative.  So,  let  us  introduce  at  this  point 
the  idea  of  a  “wavelet”: 


We  shall  speak  of  a  function  ip  in  L1(e)  as  being  normalized  and  wavelike  if 

POO  A  rjf 

(i)  ip  €  L2( M),  (u)  jf  |#0|2  y  =  1  U  €  R).  (4.2) 

By  “wavelike,”  we  are  actually  referring  to  the  characteristic  that  ip  has  average  zero; 
that  is 

POO 

I  ip(x)  dx  =  0. 

—  OO 

But  ip  6  Z4(m)  implies  that  ip  is  continuous,  so  the  existence  of  the  integral  in  (4.2)(m) 
is  enough  to  ensure  that  ip( 0)  =  0.  This  prompts  the  following  definition. 


^(0) 


(4.3)  Definition.  When  ip  is  a  normalized,  wavelike  function  in  L1( M),  the  mapping 

f  — >  T+f{r,x)  =  yy  ^  -y^)  dv  (4.3)(») 

is  called  the  continuous  wavelet  transform,  while 


is  called  the  inverse  continuous  wavelet  transform. 


dxdr 


(4.3)(») 


The  critical  idea  here  is  this:  We  devise  a  ip  with  wavelike  properties  (fulfilling 
a  role  much  like  the  complex  exponential  in  the  Gabor  transform),  but  then  we  in¬ 
tegrate  translations  and  dilations  of  this  function  against  the  signal  to  be  analyzed. 
Resultingly,  we  obtain  varying  levels  of  resolution  as  r  ranges  over  M. 

Now,  the  Gabor  transform  and  its  properties  are  very  well  known.  One  very 
important  idea  associated  with  it  is  the  Wigner-Ville  distribution.  Define  the 
“voice  transform”  Vj<g  to  be 

/OO  _ 

f(v  +  p)g(v)e2m(vq+pq/2)  dv.  (4.4) 

-OO 

By  a  simple  change  of  variables,  we  observe 

=  /”  f(v  +  dv 

J — OO 
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= rj(v+ 


(4.5) 


which  is  also  referred  to  as  the  cross  ambiguity  function  of  /  and  g.  It  is  worth 
noting  here  the  relationship  between  Vft9  and  Sgf: 

/OO  _ 

f(v  +  p)g(v)e2^v+p/2)q  dv 

-OO 

/OO  _  roo  _ 

f(v)g(v  -  p)e2m^~p^q  dv  =  e~vtpq  /  f{v)g{v  -  p)t2™q  dv 
-OO  J  —  OO 


=  e 


—mpq 


SgfiP ,  -?)• 


We  obtain  the  Wigner-Ville  distribution,  Wf#,  by  taking  the  two-dimensional 
Fourier  transform  of  V/)3  [8] : 


/CO  /«oo 

/  dqdp 

-oo  «/  —  OO 

■£ 


/(*  +  -  ^)e  27r,K 


(4.6) 


This  distribution  has  been  studied  extensively  and  has  served  as  the  prototype  of 
all  time-frequency  distributions.  [9]  Our  goal  is  to  carry  over  this  idea  to  the  wavelet 
transform.  We  take  the  following  approach: 

Represent  the  Heisenberg  group,  Ti,  associated  with  K  as  the  group  of  matrices 
of  the  following  form 


n  =  { 


For  ease  of  notation,  we  express  the  group  elements  by 


1 

X 

£ 

2s  ‘ 

> 

0 

1 

0 

1R 

0 

0 

1 

—  X 

0 

0 

0 

1 

> 

(4.7) 


1  x  £  2s 
0  1  0  ( 
0  0  1  -X 

0  0  0  1 


(4.8) 


so  matrix  multiplication  corresponds  to  the  group  operation 

(x,£,s).(y,T},t)  =  (x  +  y,Z  +  T],s  +  t  +  \{xg  -  yQ)  ■ 


(4.9) 
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It  is  quickly  verified  that  H  is  in  fact  a  group,  with  element  inverses  given  by 

(x,^,s)_1  =  -5).  (4.10) 

Now  we  introduce  the  Schrodinger  representation  tt  of  which  is  a  representation 
extending  to  H  the  basic  operations 

x(x,  0, 0)  :  f(v) — >f(v  +  x)  (x  £  IK),  (4.11) 

and 

5r(0,£,0)  :  f{v)  — »  e2*iv*f(v)  (f  €  R).  (4.12) 

These  are  the  operations  of  translation  and  modulation  which  underlie  the  Euclidean 
Fourier  transform.  Now,  since 

(x,0,0).(0,£,0).(M,0)_1.(0,£,0)_1  =  (0,0, *£)  , 

conditions  (4.11)  and  (4.12)  completely  determine  7r  on  the  subgroup  Z  =  {(0,0,  s)  : 
s  G  f }  of  7"f;  indeed,  a  simple  calculation  shows  that 

*(0,0,*fl  :/(»)—.  «**■*/(>>), 

and  so,  in  general, 

=  ?r(0,£,s  +  i^)7r(a:,0,0)/(u) 

=  tt(0,0,s  +  \x()ir(Q,(,0)f(v  +  x)  =  e^e^+^)f{v  +  *).  (4.13) 

Notice  that  this  representation  is  well-suited  to  express  the  Voice  transform,  since 

v/>,«)=  r  /(»+?)» 

J — oo 

/oo  _ 

n(p,q,0)f(v)g(v)dv.  (4.14) 

-oo 

Likewise,  the  Gabor  transform  may  be  represented  by 

Sgf(p,-q)  =  e*™Vf>g(p,q)=  f  n(p,q,^-)f(v)g(v)dv.  (4.15) 

J  —  OO  Z 

In  the  very  same  fashion,  we  represent  the  affine  group,  A ,  associated  with  1R  as  a 
group  of  matrices  of  the  following  form 

'Tlx  0  0  1 

0  r  0  0 

A  =  <  ,  ,  :x,r£IR,r>0>  .  (4.16) 

0  0  1/r  — x/r 

.[o  0  0  1  J 
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Again,  we  simplify  notation  by  expressing  the  group  elements  by 

'1x0  O' 

0  r  0  0 

(x, r) ~ 

0  0  1/r  — x/r 

0  0  0  1 


(4.17) 


The  more  common  (and  indeed  simpler)  way  to  represent  this  group  is  with  2x2 
matrices,  so  that 

,  \  n  x' 

0  r 

but  it  is  a  simple  task  to  show  that  these  are  equivalent,  and  the  form  in  Eq.  (4.17) 
is  more  appropriate  for  our  purposes  here.  Matrix  multiplication  corresponds  to  the 
group  operation 

(x,r).(y,p)  =  (px +  y,rp).  (4.18) 


The  element  inverses  in  A  are  given  by 


(x,r)  1  =  (-x/r,  1/r). 


(4.19) 


Let  us  now  incorporate  the  Heisenberg  and  affine  groups  into  a  single  representation, 
and  expand  the  nr  representation  so  as  to  include  the  wavelet  transform.  Introduce 
the  affine-Heisenberg  group,  Q ,  of  which  both  the  affine  and  Heisenberg  groups 
are  subgroups  by  letting 


'i 

X 

2s  ' 

"i 

g  =  < 

0 

0 

r 

0 

0 

1/r 

—x/r 

:i,(,s,r61,  r  >  0  > 

.0 

0 

0 

1 

> 

(4.20) 


Once  more,  it  will  be  notationally  convenient  to  write  (x,  £,  s,  r )  for  an  element  of 
Q  instead  of  its  matrix  realization.  With  this  notation,  the  group  operation  on  Q 
becomes 


(x,Z,s,r).(y,ri,t,p)  =  (px  +  y,  t/p  +  r],  s  +  t  +  \{pxy  -  y£/ p),  rp ).  (4.21) 

We  extend  the  n  representation  from  the  Heisenberg  group  to  Q  by  the  operation 

7r(0,0, 0,r)  :  f(v) — >rl!2f(rv)  (r  >  0).  (4.22) 

Consequently,  when  we  combine  Eq.  (4.22)  with  Eqs.  (4.11)  and  (4.12),  we  determine 
7r  on  all  of  Q  to  be 

ir(z,  (,  s,  r )  :  f(v)  — >  r1/2e2^e2«(^+H)y(rv  +  x)  .  (4.23) 


29 


To  verify  that  tt  is  in  fact  a  well-defined  representation  on  Q,  observe  that 

w(x,  t,  a,  r)  (ir(y,  rj,  t,p)f(v))  =  r(x,  a,  f(pv  +  y) 

—  rl/2pl/2e'*™(s+t)e2™(rvt+%xt)e2iri(p(rv+x)T)+5yv) f(p(rv  -f  x)  +  l /). 

On  the  other  hand, 

*((*,£, s,r).(y,r),t,pj)f(v) 

=  tt (px  +  y,  (/p  +  i 7,  s  +  t  +  \(pxr)  -  y£/p),  rp)f(v ) 

—  ^rpy/2e2'ni(s+t+l(P^V-yi/p)) e2iri(rpv(i/p+ri)  +  ^(px+y)(C/p+j,))  ^  _|_  pX  _|_  y) 

=  r1/2^1/2e2iri(s+<)e27rt'(r^+|-x^)e27ri(p(rt/+i)))+|OT) f(p(rv  +  x)  -)-  y). 

Hence,  n(x,  £,  s,  r)  q ,  t ,  p)/(u))  =  tt  ((x,  £,  s,  r).(y,  p,  i,  p))f(v). 

Now  we  finally  have  a  representation  which  puts  the  Gabor  and  wavelet  trans¬ 
forms  “under  one  roof,”  so  to  speak.  Observe  that  in  the  case  s  =  ^  and  r  —  1 ,  Q 
reduces  to  the  Heisenberg  group  and  by  Eq.  (4.15)  the  Gabor  transform  is  represented 
as 

SgfiP,  ~q)  =  [  *(p,  q,  i?,  0 )f(v)g(v)  dv. 

J — oo  Z 

On  the  other  hand,  letting  £  =  0  and  s  =  0,  we  obtain 

/CO  _  POO  _ 

7r(x,  0, 0,  r)f(v)^(v)  dv  =  r1'2  /  f(rv  +  v)  dv 

-oo  J  — oo 

=  ^  =  Tii>f{r,x), 

which  gives  us  a  representation  of  the  wavelet  transform. 

So,  let  us  refer  back  to  the  Wigner-Ville  distribution  and  try  to  find  some  analo¬ 
gous  setting  in  the  wavelet  domain.  Recall  that 

wiAxi  0  =  f{vfAPiq))  =  *?r(^oo7r(p)?,o,o)/(u)(7(u)d^, 

where  T  here  denotes  the  two-dimensional  Fourier  transform  on  1R2.  In  this  case,  Vfi9 
takes  signal  /  and  window  g  to  a  function  on  the  plane  {(p,  q)  :  p,q  E  M}.  The  group 
structure  with  regard  to  both  variables  p  and  q  is  additive ,  hence  T  is  the  “natural” 
transform  to  apply  to  Vj>g. 

Now  consider 

/OO  _ 

7r(x,0,0  ,r)f(v)ip(v)dv. 

-oo 
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Here  T^f  takes  signal  /  and  wavelet  %f>  to  a  function  on  the  upper  half-plane  {(r,  x)  : 
r,  x  €  M,r  >  0}.  But  we  observe  that  r  plays  a  scaling  role  here,  and  hence  the  group 
structure  with  regard  to  this  variable  is  multiplicative.  As  we  have  seen,  the  action  of 
the  Mellin  transform  on  the  multiplicative  group  is  analogous  to  that  of  the  Fourier 
transform  on  the  additive  group.  This  leads  us  to  consider  the  effects  of  the  Mellin 
transform  on  the  scale  variable  in  T^f  (we  will  utilize  both  and  A4+  here): 


Let  /  and  xjj  be  in  H\{ R).  Then 


M:T^f — >  jf  ^ I  7r(x,0,0 ,r)f(v)^(v)dv^r  2lTt 
=  J  ^ j  rx^2f{rv  -f  x)^>{v)  dv^jr~2mt~l  dr. 


By  Parseval’s  theorem,  and  since  /,  6  7/+(R), 


I  (y  rl/2f(rv  +  x)ip(v)  dvjr  2mt  1 dr 

= ?■ 

Now  we  change  the  order  of  integration  and  make  the  variable  change  u  =  £/r  to 


obtain 


-2-Kit  dr 


if  QT  r-V2e2^/r)^0^^  d^r~27rit  y 

=  (jT°  S0r2^<_1/2  (^°°  e2*ixuf(u)u2lTit-l/2  du 

=  K${t)(Tx*Kf){t) 


(since  /C  =  .M4.  o  T  for  functions  in  77?  (m)).  Thus, 


:  T/,/  — >  (JCil>)(Tx  *  ICf). 


(4.24) 


This  result  is  of  particular  interest  for  at  least  two  reasons.  First,  notice  that 
taking  the  Mellin  transform  with  respect  to  the  scale  variable  in  T^f  allows  us  to 
“factor  out”  the  wavelet  transform,  i.e.,  9 Jl(T^f)  is  a  product  of  a  transform  of  xf>  with 
a  transform  of  /.  Second,  Eq.  (4.24)  allows  a  representation  of  T^,f  which  is  entirely 
analogous  to  the  standard  representation  of  the  Wigner-Ville  distribution.  Indeed, 
taking  the  Fourier  transform  of  the  voice  transform  Vjt9  with  respect  to  both  variables 
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we  obtained  the  Wigner-Ville  distribution,  most  commonly  written,  as  in  Eq.  (4.6), 
as 

W,Jx,0  =  jT  /(*  +  rfp.  (4.25) 

Analogously,  taking  the  Mellin  transform  with  respect  to  the  scale  variable  and  the 
Fourier  transform  with  respect  to  the  shift  variable  of  the  wavelet  transform  T^f  we 
obtain 

/OO 

m{T4,f)(t,x)e-2™*dx 

-OO 

/OO  _ 

(%*JCf){t)K^{t)e-2mxidx 

-oo 

Because  this  mapping  is  exactly  the  affine  analogue  of  the  Wigner-Ville  distribution, 
we  shall  call  the  function 

WfAW)  =  (/ (4.26) 

the  affine  Wigner-Ville  distribution  of  /  and 

Previous  attempts  at  affine  versions  of  the  Wigner-Ville  distribution  have  been 
made  by  Altes[10],  Flandrin  [11],  and  Parks  and  Shenoy  [8],  but  none  utilize  the 
V- transform.  What  we  have  done  is  begin  an  analysis  of  the  ^-transform  in  a  setting 
which  unifies  both  the  Gabor  transform  and  the  wavelet  transform.  Its  application  to 
this  unified  setting,  and  its  suitability  to  the  work  of  Altes  mentioned  in  the  previous 
section,  indicate  that  the  /C-transform,  as  the  composition  of  Fourier  and  Mellin 
transforms,  is  perhaps  a  worthwhile  topic  for  further  investigation.  In  fact,  it  appears 
this  transformation  may  be  an  underlying  element  in  the  realm  of  wavelets  and  other 
such  methods  of  time-scale  analysis;  as  these  areas  are  further  pursued  so  might,  also 
the  /C- transform  be  more  thoroughly  considered. 
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